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Abstract 

Molecular dynamics simulations are carried out to investigate the dynamic behavior of the slip 
length in thin polymer films confined between atomically smooth thermal surfaces. For weak wall- 
fluid interactions, the shear rate dependence of the slip length acquires a distinct local minimum 
followed by a rapid growth at higher shear rates. With increasing fluid density, the position of the 
local minimum is shifted to lower shear rates. We found that the ratio of the shear viscosity to the 
slip length, which defines the friction coefficient at the liquid/solid interface, undergoes a transition 
from a nearly constant value to the power law decay as a function of the slip velocity. In a wide 
range of shear rates and fluid densities, the friction coefficient is determined by the product of the 
value of surface induced peak in the structure factor and the contact density of the first fluid layer 
near the solid wall. 
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I. INTRODUCTION 



The fluid flow in microfluidic channels can be significantly influenced by liquid slip at the 
solid boundary because of the large surface to volume ratio The Navier model for the 
partial slip boundary conditions relates the fluid slip velocity to the tangential viscous stress 
at the wall via the friction coefficient, which is defined as the ratio of the shear viscosity to 
the slip length. Geometrically, the slip length is determined as a distance from the solid wall 
where the linearly extrapolated fluid velocity profile vanishes. Experimental studies have 
demonstrated that for atomically smooth surfaces the slip_^length is relatively large for non- 



wetting liquid/solid interfaces [2,0,141, polymeric fluids j5|,l6|] and high shear rates 



,y,3. 



By 



contrast, surface roughness even on the submicron length scale strongly reduces the degree 
of slip for both polymeric [lo| and Newtonian 11, 1^. If^ fluids. However, experimentally 
it is often difficult to isolate and, consequently, to analyze the effects of wettability, surface 
roughness and shear rate on slip boundary conditions because of the small scales involved. 
■ 3uring the past two decades, a number of molecular dynamics (MD) simulations 
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23| have been performed to investigate the dependence of the slip 



length on the structure of dense fluids near atomically flat surfaces. For simple fluids, it has 
been well established that the slip length correlates inversely with the wall-fluid interaction 



energy and the surface induced order in the first fluid layer 
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23| . The strong wall- 



fluid coupling and commensurate structures o: 



boundary condition for monoatomic fluids 15 



;he liquid/solid interface lead to the no-slip 



2l| . The interfacial slip is enhanced in the 



flow of polymer me 



the flat wall 
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ts due to the higher shear viscosity and reduced molecular ordering near 



I, |26| . The magnitude of the slip length in the flow of polycarbonate 
blends past atomically smooth nickel surfaces is determined by the competition between 
the apparent slip of the adsorbed polymer layer and its thickness 27(|. The slip boundary 
conditions for the polymer melt flow over the brush can be controlled by changing the density 
of the end-grafted polymer layer Finally, the degree of slip at the interface between 

immiscible polymers is larger for longer chain molecules and greater repulsion between unlike 



species 



29|. 



Until ten years ago, the slip length in the Navier model was assumed to be independent 

I P 

of shear rate. In the seminal MD study by Thompson and Troian [30] on shear flow of simple 
fluids past atomically flat rigid walls, it was found that the slip length increases nonlinearly 
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with the shear rate. Upon increasing the surface energy, a gradual transition from the 



power law behavior [30] to the linear rate dependence was recently reported 
fluid interactions and incommensurate structures of the liquid/solid interface 



or weak wall- 



231] . In a wide 



range of shear rates and surface energies, the slip length in simple fluids is determined by a 
function of the in-plane structure factor, contact density and temperature of the first fluid 
layer near the wall [23I. The rate-dependent slip length at the interface between a polymer 
melt composed of short linear chains and a crystalline substrate is also well described by the 
power law function used to fit the data for simple fluids 0, 3^. Whether these boundary 
conditions hold for higher molecular weight polymers or rough surfaces is an open question 
that has to be addressed for the correct interpretation of the experimental results and 
modeling fluid flows in microfluidic channels. 

In this paper, the effects of shear rate and fluid density on slip boundary conditions for 
the flow of a polymer melt past a crystalline wall are investigated by molecular dynamics 
simulations. We will show that the rate-dependent slip length reaches a local minimum 
and then increases rapidly at higher shear rates. In the shear thinning regime, the friction 
coefficient at the liquid/solid interface follows a power law decay as a function of the slip 
velocity. In a wide range of shear rates and polymer densities, the ratio of the slip length 
to the shear viscosity is determined by the surface induced structure in the first fiuid layer. 
These results suggest that the correlation between the fiuid structure near the fiat wall and 
the slip length is determined by a universal combination of parameters for both polymer 
melts and simple fiuids [23|]. 

This paper is organized as follows. The molecular dynamics simulation model and pa- 
rameter values are described in the next section. The results for the rate dependence of 
the slip length and the analysis of the fiuid structure near the solid wall are presented in 
Section IIIII The conclusions are given in the last section. 



II. DETAILS OF MOLECULAR DYNAMICS SIMULATIONS 



The polymer melt is confined between two atomically fiat walls and is subject to planar 
shear in the x direction (see Fig. [1] for a schematic representation). The fiuid monomers 
interact through the pairwise Lennard- Jones (LJ) potential 

. r / V r . 



VLj{r)=Ae 



(1) 



with a cutoff distance Vc = 2.5 a, where e and cr represent the energy and length scales of 
the fluid phase. The total number of monomers is Nf = 6000. The wall-fluid interaction 
parameters of the L J potential are set to e„{ = 0.9 e and cr^f = throughout the study. 
The polymer melt is modeled as a collection of bead-spring chains of = 20 L J monomers 
interacting through the FENE (finitely extensible nonlinear elastic) potential 31 1 

VrENE{r) = -'^r^Hl-ryr^], (2) 

with Kremer and Grest parameters = 30 ea^'^ and Tq = 1.5 cr [32]. This choice of parameters 



prevents polymer chains from unphysical crossing or bond breaking c 



33\ 



The system was coupled to an external heat bath through a Langevin thermostat 
with a random force and a damping term with a small value of the friction coefficient 
r = 1.0r~^ [34^. The thermostat was only applied to the equation of motion for a fluid 
monomer along the y axis to avoid a bias in the shear flow direction [l^. The equations of 
motion for the x, y and z components are given by 

dxi 



= - J]— ^, (3) 



Qy. . 

mjji + mVyi = - + fi , (4) 

where fi is a random force with zero mean and variance (/i(0)/j(t)) = 2mkBTT5{t)5ij, the 
temperature is T =l.le/kB and is the Boltzmann constant. The equations of motion 
for fluid monomers and wall atoms were integrated using the fifth-order gear-predictor al- 
gorithm |35] with a time step At = 0.002 r, where r= ^Jma'^ /e is the characteristic time of 
the LJ potential. 

Each wall was constructed out of 576 LJ atoms, which formed two (111) layers of the 
face-centered cubic (fee) lattice with the density = 1.40 a^'^. The corresponding nearest- 
neighbor distance between equilibrium positions of the wall atoms in the xy plane is c/ = 1.0 cr. 
The dimensions of the Couette cell {xyz) were held fixed to 20.86 cr x 12.04 a x h, where h is 
the channel height (see Table[T]). The MD simulations were performed at a constant density 
ensemble. The fluid density is defined as a ratio of the total number of monomers to the 
volume accessible to the fluid, i.e., 0.5 cr„ away from the fee lattice planes in the z direction. 
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The range of fluid densities considered in this study is 0.86 ^ pcr^ ^ 1.02. The pressures of 
the system at equihbrium are given in Table[T] for each value of the fluid density. Periodic 
boundary conditions were employed along the x and y directions parallel to the confining 
walls. 

The wall atoms were attached to their equilibrium fee lattice positions by the harmonic 
spring with the potential Vsp = \ nr"^. The spring stiffness k = 1200 e/cr^ is large enough so 
that the mean-square displacement of the wall atoms is less than the Lindemann criterion 
for melting, i.e^{5v?) / d"^ < 0.023 (see Ref. |36j). In the previous paper on shear flow of 
simple fluids [37|, it was shown that the slope of the rate-dependent slip length remains 
approximately the same for the thermal walls with k = 1200 e/cr^ and rigid fee walls. The 
friction term and the random force were also added to the x, y and z components of the 
wall atoms equations of motion 

Xi + m^Txi = - ^ - + fi , (6) 

m + mu,Tyi = - ^ - + fi , (7) 

Zi + mujTzi = - ^ - + fi , (8) 

where the mass of the wall atoms is m^^, = 10 m and the summation is performed only over 
the fluid monomers within the cutoff distance rc = 2.5 a. 

The simulations began with an equilibration period of about 5 x lO'^r with both walls being 
at rest. Then, the steady-state shear flow was generated by translating the upper wall with a 
constant velocity U in the x direction, while the lower wall always remained stationary (see 
Fig.[l]). The fluid velocity proflles were computed by averaging the instantaneous monomer 
velocities in the x direction within bins of thickness Az = 0.2 a for a time interval up to 

4 X lO^r at the lowest shear rates. The density proflles were averaged for a time period 

5 X lO^'r within bins of thickness of only Az = 0.01 a to accurately resolve the fluid structure 



near the walls 



231]. The maximum value of the Reynolds number is i?e ~ 10.5, which was 



estimated from the maximum fluid velocity difference across the channel, the shear viscosity, 
and the channel height. The small values of the Reynolds number correspond to laminar 
flow conditions even at the highest shear rates considered in this study. 
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III. RESULTS 



A. Fluid density and velocity profiles 

The fluid structure in the direction perpendicular to a flat surface usually consists of 



several distinct layers on the molecular scale [38|. Examples of the averaged monomer 
density profiles near the lower thermal wall are presented in Fig. [2] for three values of the 
upper wall velocity. The highest peak is associated with a layer of fluid monomers (each 
connected to a polymer chain) in contact with the wall atoms. The contact density Pc 
corresponds to the maximum value of the peak. The density oscillations gradually decay 
away from the wall to a uniform profile with bulk density p = 0.91(T^^. The amplitude of 
the density oscillations is reduced at higher values of the slip velocity. 

FigureO shows the averaged velocity profiles for the same values of the upper wall speed 
as in Fig. [2] and the fiuid density p = 0.91(T^'^. For a given U, the relative slip velocity is 
the same at the upper and lower walls. Surprisingly, the scaled slip velocity is smaller at 
the intermediate speed of the upper wall U = 1.0 a/r. Because of the thermal fiuctuations, 
longer averaging times (up to 4 x lO^r) were required to accurately resolve the fiuid velocity 
profile for U = 0.05 a/r. A small curvature in the velocity profile is evident for the largest 
velocity of the upper wall U = 5.0 a /t. The deviation from the linearity might be related to 
the heating up of the fiuid near the walls (see discussion below). The slope of the velocity 
profile (used for calculation of the slip length and the shear rate) was computed from the 
central part of the velocity profile excluding the data from the region of about 2 a near the 
walls. 

We comment that at small shear rates, 'jr < 0.01, and higher fiuid densities, p > 1.04 a~^, 
the fiuid velocity profiles acquired a pronounced curvature within a distance of about A a 
from the walls. The linearity of the fiuid velocity profiles is restored at higher shear rates. 
This effect will be examined in a separate study. 

B. Rate dependence of the slip length and shear viscosity 



In the MD study by Thompson and Troian 301] on shear fiow of Newtonian liquids past 
atomically fiat surfaces, the slip length was found to increase nonlinearly with the shear 
rate. The data for different wall densities and weak wall-fluid interactions were fltted well 
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by the power law function 

^.(7) = ^:(l-7/7c)-°^ (9) 

where 7c and L° are the parameters used for normahzation. The rate-dependent shp length 
in the Poiseuille flow of simple fluids was also well described by the function [Eg. (El)] for 



incommensurate structures of the liquid/solid interface and weak surface energy 23j. 

The nonlinear behavior of the slip length can be interpreted in terms of the velocity- 
dependent shear stress at the wall [391]. By definition, the fluid slip velocity K is related to 
the slip length as K = "^Lg. On the other hand, the shear stress through any plane parallel 
to the walls is the same in the steady-state shear flow, and it is equal to the fluid viscosity 
multiphed by the shear rate, i.e., axz = l^i- At the liquid/solid interface the viscous shear 
stress becomes a^z = kVg, where k is the friction coefficient per unit area. The functional 
form given by Eq. Q) implies that the friction coefficient, k= ^/Ls, can be expressed as a 
function of the slip velocity as follows: 

/c(K) = Ci(v/C2 + K'-K), (10) 



where Ci = /i/27c(L°)^ and Co = (2L°7c)^, and /i is the shear-independent viscosity [20|, l30|. 



In the previous MD study 20|], it was shown that the slip length increases monotonically 
with shear rate for polymer melts with short linear chains ^ 16. The data for the slip 
length were well fitted by the power law function given by Eq. ([9]) for shear rates exceeding 
7 > 5 X 10^'^ r^^. Although the fluid viscosity exhibited a slight shear thinning behavior for 
chain lengths ^ 8, no minimum in the rate dependence of the slip length was observed. 
In the present study, the shear thinning effect is more pronounced because the polymer melt 
consists of longer chains and the simulations are performed at higher fluid densities. 

FigureHlshows the rate behavior of the fluid viscosity, which was computed from the Kirk- 
wood relation for the shear stress [si!] in steady-state shear flow. In the range of accessible 
shear rates, the transition from the Newtonian regime to the shear thinning is more evident 
at lower fluid densities. At higher shear rates, the bulk viscosity approximately follows a 
power law decay with the exponent —0.37 (see Fig.Sj). A similar shear thinning behavior of 



4Q|, 



4l|. The error bars are larger 



polymer melts was reported in the previous MD studies 
at lower shear rates because of an increasing role of thermal fluctuations 421]. 

The dynamic response of the slip length with increasing shear rate is shown in Fig. [5] for 
the indicated values of the polymer density. The rate dependence of the slip length exhibits 



a distinct local minimum, which is followed by a rapid increase at higher shear rates. The 
minimum is shifted to lower shear rates and becomes more shallow at higher fluid densities. 
The slip length tends to saturate to a constant value at the lowest shear rates for the fluid 
density p = 0.86a^^. Typical error bars are included for three points corresponding to the 
fluid velocity profiles shown in Fig. [31 In the range of shear rates from the local minimum 
up to the highest 7, the data for the slip length (as well as the ratio Ls/ p) cannot be fitted 
by the power law function given by Eq. ([9]). For each value of the fluid density, the highest 
shear rate corresponds to the maximum sustainable shear stress at the interface 23 1. 

The data for the rate-dependent slip length and fluid viscosity are replotted in Fig. [6] 
in terms of the friction coefficient at the liquid/solid interface and the slip velocity. The 
constant value of the friction coefficient at small slip velocities indicates that the ratio of 
the slip length and viscosity is rate-independent. At higher fluid densities (or pressures), 
the friction coefficient is larger, and the transition to the power law regime is shifted to 
smaller slip velocities. The transition point approximately corresponds to the location of 
the minimum in the shear rate dependence of the slip length shown in Fig.O The data for 
the lowest fluid density p = 0.86 a"^ can be fairly well fitted by Eq. ( ITOl) [see Fig. [6]], but the 
fit becomes worse at higher fluid densities (not shown). The simulation results indicate that 
the nonmonotonic rate behavior of the slip length for polymer melts can be ascribed to the 
velocity dependence of the friction coefficient, which undergoes a transition from a nearly 
constant value to a decreasing function of the slip velocity. 

Furthermore, the data for different fluid densities can be collapsed onto a single master 
curve by scaling both axes (see Fig.[71). The values of the normalization parameters for 
the friction coefficient and the slip velocity are listed in the inset. The decay at large slip 
velocities approximately follows the power law function with the exponent —0.7 (see Fig.[7j). 
The onset of the power law decay depends on the polymer density and might occur when 
the time scale of the in-plane diffusion of the fluid monomers over the lattice spacing d is 
comparable to d/Vg. We also note that the data for the viscosity and slip length at the 
interface between a polymer melt (with linear chains = 20 and density p = 0.88 o"^^) and a 
rigid fee wall (with density = 1.94 and £„f = e) can be normalized to follow the master 
curve shown in Fig. [7] 43 1. 
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C. Analysis of the fluid structure in the first layer 



The surface induced order in the adjacent fluid 
at the reciprocal lattice vectors of the crystal wall 



ayer is described by the structure factor 



15|. The structure factor is defined as 



A k r. 



(11) 



where = {xj, yj) is the position vector of the j-th monomer and the summation is performed 
over Ni monomers within the first layer. The friction force between a layer of adsorbed 
monomers and a solid substrate is proportional to the peak value of the in-plane structure 
factor estimated at the first reciprocal lattice vector 4J, |45|] . Previous MD studies of simple 
fluids have also reported a strong correlation between the slip length (which is inversely 
proportional to the friction coefficient at the liquid/solid interface) and the surface induced 



order in the first fiuid layer 



3,[l6, 



23| 



Figured] shows the averaged structure factor in the first fiuid layer moving with a constant 
velocity against the stationary lower wall. The circular ridge at the wavevector |k| ^271/ a 
is attributed to the short-range ordering of the fiuid monomers. The periodic potential of 
the crystalline substrate induces several sharp peaks in the fiuid structure factor, which 
are reduced at higher slip velocities. The first reciprocal lattice vector in the direction of 
shear fiow is Gi = (7.23 cr^^, 0). Although the nearest-neighbor distance between equilibrium 
lattice sites in the xy plane is equal to the diameter of the LJ monomers, the surface induced 
order in the first fiuid layer is frustrated by the topological constraints associated with 
packing of polymer chains near the surface. 

In the recent MD study on shear fiow of simple fiuids past rigid fee walls [23], it was 
shown that the slip length scales according to 



:Ti/5(Gi)p, 



1.44 



(12) 



where Gi is the the first reciprocal lattice vector, Ti and pc are the temperature and the 
contact density of the first fiuid layer respectively. This relation was found to hold in a wide 
range of shear rates and for weak wall-fiuid interactions 0.3 ^ ^wf/ £ ^ l-l- It is expected, 
however, that the scaling relation given by Eq. (1121) is not valid at higher surface energies 
because the first fiuid layer is locked to the solid wall and, as a consequence, the slip length 
is „e,a«ve No.e ..t t.e surface .duced peaUs s..uc..e faeto. are much higher 



for polymers (see Fig. [8]) than for simple fluid [23|. Since the viscosity is independent of 
shear rate for simple fluids, the scaling relation Eq. f|T2l) also predicts the dependence of the 
friction coefficient on microscopic parameters of the adjacent fluid layer. 

We computed the parameters entering Eq. f[T^ in the first fiuid layer for the polymer 
melt. The fiuid temperature was estimated from the kinetic energy 

N 

3N 



kBT=^Y.ir,-Hr^)]^ (13) 



i=l 

where is the instantaneous velocity of the fiuid monomer and v(rj) is the local average 
fiow velocity. Averaged temperature profiles are shown in Fig. |9] for the same fiow conditions 
as in Figs. [2] and [3l At low shear rates 7 < 0.01 r~^, the fiuid temperature is equal to its 
equilibrium value ofT = l.le/kB determined by the Langevin thermostat. As in the case of 
simple fiuids {23! . Is?! , with further increase of the shear rate, the polymer melt heats up and 
the temperature profile becomes non-uniform across the channel. The fiuid temperature is 
higher near the walls because of the large slip velocity, which is comparable to the thermal 
velocity, v^= ksT/m, at high shear rates. The non- uniform temperature distribution might 
be related to the slight curvature in the velocity profile for the top wall speed U = 5.0 a /t 
shown in Fig.O We can also notice that at high shear rates, the temperature in the y 
direction, in which the Langevin thermostat is applied, is slightly smaller than its value in 
the X and z directions (see inset of Fig. [9]). This difference implies that the kinetic energy 
in the y direction dissipates faster than the energy transfer from the other directions. In 
principle, it is important to investigate whether the choice of the thermostat, e.g. dissipative 
particle dynamics thermostat, has an effect on the velocity and temperature profiles at high 
shear rates [46|. In the present study, the temperature in the first fiuid layer was computed 



from 

Ti= r T{z)p{z)dz / r p{z)dz, (14) 

J ZQ J Zo 

where the limits of integration {zq = — 14.4cr and zi = — 13.8cr) are given by the width of 
the first peak in the density profile (e.g. see Fig. [2]). 

The ratio of the slip length to viscosity is plotted in Fig.[Tn]as a function of Ti/[S{Gi) pj 
for different fluid densities and shear rates. Except for higher densities at low shear rates, 
the data collapse onto a single curve, which approximately follows the power law function 
with a slope of 0.9 (see Fig. fTOl) . Although the slope is different from the value 1.44 reported 
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for simple fluids 23(], it is remarkable that the same combination of parameters determines 
the boundary conditions for both polymer melts and simple fluids. The simulation results 
also show that for each value of the fluid density, the number of fluid monomers in the flrst 
layer Ni = S{0) decreases by about 4% at the highest shear rates. The same data for the slip 
length and viscosity are replotted as a function of the variable 5'(0)/ [5'(Gi) pd in Fig.fTTl 
The collapse of the data is consistent with the previous results for the rate-independent 
slip length at the interface between simple fluids and crystalline walls The correlation 
between the friction coefficient and the fluid structure in the adjacent layer presented in 
Figs. [To] and fTTl suggests a possible existence of a general functional relation between the slip 
length, shear rate and polymer density. 



IV. CONCLUSIONS 



In this paper, the rate dependence of the slip length in the shear flow of polymer melts 
past atomically smooth, thermal surfaces was studied using molecular dynamics simulations. 
For weak wall-fluid interactions, the slip length passes through a minimum as a function of 
shear rate and then increases rapidly at higher shear rates. The nonlinear rate dependence 
of the slip length was analyzed in terms of the dynamic friction at the liquid/solid interface. 
In a wide range of fluid densities, the friction coefficient (the ratio of the shear viscosity 
and the slip length) undergoes a universal transition from a constant value to the power law 
decay as a function of the slip velocity. The simulation results for polymer melts confirm the 
previous findings for simple fiuids |23|] that the surface induced structure and the contact 
density of the adjacent fluid layer are crucial factors for determining the value of the friction 
coefficient. 

Future research will be focused on the rate dependence of the slip length for higher fluid 
densities and entangled polymer melts. 
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FIG. 1: (Color online) A schematic representation of the steady-state shear flow in the Couette 
cell. The upper wall is translated with a constant velocity U in the x direction. The fluid slip 
velocity is determined from the relation Vs = jLg, where 7 is the slope of the velocity profile and 
Lg is the slip length. 
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FIG. 2: (Color online) Averaged fluid density profiles near the stationary thermal wall with 
K = 1200e/a^ and e„f/e = 0.9. The velocities of the upper wall are tabulated in the inset. The 
fluid density in the center of the cell is p = 0.91 cr~3. The left vertical axis denotes the position of 
the liquid/solid interface at 2:= — 14.74 ex. 
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FIG. 3: (Color online) Average normalized velocity profiles, Vx/U, for the indicated values of 
the upper wall velocity U and the fluid density p = 0.91(T~^ . Vertical axes denote the position of 
the fee lattice planes at z= — 15.24 a and z = 12.03 o". The dashed lines indicate the location of 
liquid/solid interfaces at a distance 0.5 (T„ away from the fee lattice planes, i.e., z = — 14.74(7 and 
z = 11.53 fj. The inset shows the same data near the lower wall. 
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FIG. 4: (Color online) Behavior of the fluid viscosity /i/ era^'^ as a function of shear rate for the 
tabulated values of the polymer density. The straight dashed line with a slope —0.37 is shown for 
reference. Solid curves are a guide for the eye. 
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FIG. 6: (Color online) Log-log plot of the friction coefficient, k= fi/Lg, as a function of the slip 
velocity for the indicated values of the fluid density. The dashed curve is the best fit to Eq. (jlOp 
with Ci = 1.77 mo""^ and C2 = 0.71 ct^t~^. Solid curves are a guide for the eye. 
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FIG. 7: (Color online) Master curve for the friction coefficient as a function of the slip velocity. The 
same data as in Fig.[6l The values of the friction coefficient k* [ era^'^] and the slip velocity V^* [ c/r ] 



used for normalization of the data are shown in the inset. The dashed line y = [1 + (4x 
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IS 



plotted for reference. 
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FIG. 8: Structure factor S{kx,ky) in the first fluid layer for the fluid density /3 = 0.91cr ^. The 
upper wall velocity \sU = O.OSa/r (a) and U = B.Oa/r (b). The same flow conditions as in Fig.O 
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FIG. 9: (Color online) Temperature profiles for the indicated values of the upper wall speed U and 
the fluid density p = 0.91 o"^^. Vertical dashed lines denote the position of the liquid/solid interface. 
The inset shows the x, y and z components of the temperature profile near the stationary lower 
wall when the velocity of the upper wall is U = S.Ocr/r. 
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FIG. 10: (Color online) Log-log plot of Ls//x as a function of the ratio Tiks/e [S{Gi) PcO"^] ^ for 
the indicated values of the polymer density. The solid line with a slope 0.9 is plotted for reference. 
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FIG. 11: (Color online) Log- log plot of the ratio Ls/fJ, as a function of the variable 
S (0) / [S (Gi) pccr^]- The values of the fluid density are tabulated in the inset. The same data 
as in Fig.[TOj The straight solid line with a slope 1.15 is plotted as a reference. 
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TABLE I: The fluid pressure P at equilibrium, i.e., with the stationary upper waU, and the channel 
height h (defined as a distance between the inner fee planes) as a function of the fluid density. 
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